
AIAA 94-1579 

Thermoelastic Formulation of Stiffened, 
Unsymmetric Composite Panels for Finite 
Element Analysis of High Speed Aircraft 


Craig S. Collier 

Lockheed Engineering and Sciences Company 
NASA Langley Research Center 
Hampton, Virginia 

Work performed under NASA contract 


AIAA/ASME/ASCE/AHS/ACS 
35th Structures, Dynamics, & 
Materials Conference 

April 18-20, 1994/Hilton Head, SC 

For permission to copy or republish, contact the American Institute of Aeronautics and Astronautics 

370 L’Enfant Promenade, SW • Washington, D.C. 20024 



Thermoelastic Formulation of Stiffened, Unsymmetric Composite Panels for 
Finite Element Analysis of High Speed Aircraft 


Craig S. Collier* 

Principal Structures Engineer 
Lockheed Engineering and Sciences Co. 
NASA Langley Research Center 
Hampton, VA 


ABSTRACT 


An emerging technology need for capturing 3-D panel 
thermoelastic response with 2-D planar finite element 
models (FEMs) is aided with an equivalent plate stiffness 
and thermal coefficient formulation. The formulation is 
general and applies to all panel concepts. Included with 
the formulation is the ability to provide membrane-bending 
coupling of unsymmetric sections and calculation of all 
thermal expansion and bending responses from in-plane 
and through-the -thickness temperature gradients. Thermal 
residual strains for both the laminates and plies are 
included. The general formulation is defined and then 
applied to a hat-shaped, corrugated stiffened panel. 
Additional formulations are presented where required to 
include all of the hat’s unique characteristics. Each 
formulation is validated independently with 3-D FEA. 


NOMENCLATURE 

laminate or panel 

AT In-plane temperature gradient 

AG Through-the-thickness 

temperature gradient 

(i = l 2, 3) 

h Distance from the reference plane 

Ajj, Bjj, Membrane, membrane-bending coupling, 

D and bending stiffness terms 

Aj a , Bi a , Membrane, membrane-bending coupling, 

D,“ and bending thermal force & moment 

coefficients 

Qjj Reduced transformed laminate elasticities 

Reduced transformed laminate thermal 
force coefficients 


* Senior member AIAA 

This paper is a work of the U.S. Government and is not subject to 
copyright protection in the United States. 


(i = x, y, xy) 

Nj, Mj Reference plane forces and moments 

e ; , Kj Reference plane strains and curvatures 


superscripts 

P 

T 

M 

m 


Panel 
Thermal 
Mechanical 
exponent value 


INTRODUCTION 


Multidisciplinary conceptual design of high speed aircraft 
requires a quick structural analysis capability for vehicle 
optimization. High speed aircraft are often designed with 
stiffened and sandwich panels fabricated from fiber- 
reinforced composite materials. The need for a quick 
analysis and the complexity of composite material, 
stiffened panels encourages approximations in the formula- 
tion of panel stiffness terms. The accuracy of approximate 
panel formulation is further diminished for hot environ- 
ments caused by supersonic flight. Temperature gradients 
induce forces and moments which must be quantified with 
thermal expansion and bending coefficients. Formulation 
of thermal coefficients is also complex, so they too, are 
usually approximated. 

This paper describes an equivalent plate stiffness and 
thermal coefficient formulation that applies to any panel 
concept, see Fig. 1. All stiffness and thermal coefficient 
data for membrane [A] & {A“}, bending [D] & {D“}, and 
unsymmetric membrane-bending coupling [B] & {B“} are 
captured for any stiffened shape. This paper examines the 
formulation's ability to analyze for the most general 
mechanical and thermal loading conditions by comparing 
its results with those produced by discrete, three dimen- 
sional FEA. The applicability of plate and beam analogies 
are addressed as well as the Kirchoff hypothesis, plane - 
sections-remain-plane. 
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Corrugated-shaped, 
stiffened panel 


Hat Stiffened 





Waffle Grid 



J Stiffened 



Honeycomb Sandwich 



Trusscore Sandwich 


Fig. 1 The formulation can be applied to any stiffened 
or sandwich, composite panel concept. 


Stiffened Panels 

Stiffened and sandwich panels efficiently provide buckling 
stability. They can carry more service load than 
unstiffened plates for a given unit weight. Stiffened panels 
are quite efficient for lightly loaded areas and applications 
of high temperature gradients. These qualities make them 
desirable for use as hot structure on high speed vehicles 
where weight reduction is a paramount objective. 

Stiffened panels are unsymmetric, primarily by nature of 
their shapes as seen if Fig. 1. Unsymmetric stiffened 
shapes produce coupling even when the panel is fabricated 
with conventional isotropic materials. Sandwich panels 
become unsymmetric due to the effect of through-the- 
thickness temperature gradients on their material elastici- 
ties. Quantification of unsymmetric behavior is important 
because it significantly alters panel response. 
Unsymmetric behavior causes coupling between membrane 
and bending panel response. Therefore, bending will 
either shorten or lengthen the panel midplane. Likewise, 
a change in panel length will create curvature. This 
membrane-bending coupling is quantified in classical 
lamination theory [1] with the [B] stiffness matrix. If 
change in panel shape and size is due to a change in 
temperature, then corresponding thermal coupling 
coefficients noted as {B“} must also be quantified. 

A method for accurately including composite lamina and 
laminate data in the formulation of stiffened panel 
structural properties was first presented in reference 2 and 


later published in reference 3. Thermal coefficients created 
to handle both in-plane and through-the-thickness tempera- 
ture gradients for membrane, bending, and membrane- 
bending coupling were also introduced in reference 2 and 
then mathematically proven in reference 3. Another 
paper, reference 4 showed how to input these data into the 
MSC/NASTRAN™ finite element analysis (FEA) program 
[5] using a model with a single plane of finite elements. 
The significance of including them for an entire aircraft 
analyses is reported in reference 6 and 7. 

Finite Element Analysis 

In-plane and through-the-thickness temperature gradients 
may be correctly applied and solutions obtained for 
anisotropic/orthotropic, unsymmetric, and unbalanced 
laminates or stiffened panels using a single plane of shell 
elements with the MSC/NASTRAN FEA program [3,4], 
This is accomplished by including the full complement of 
smeared equivalent plate stiffness matrices and thermal 
expansion and bending coefficient vectors in the FEM. 

Smeared equivalent plate stiffness and thermal coefficient 
formulations of this paper are particularly useful for 
coarsely meshed models of a complete structural entity 
such as an engine or airframe. Practical models of such 
large surface areas can only be accomplished with a single 
plane of shell finite elements. Too many elements are 
necessary to construct a discrete three-dimensional panel 
model over a complete airframe surface. 



Missing: 

• A12 & D12 directional coupling terms 

• (1-v 2 ) plate poissons term 

• Unsymmetric data 



Missing: 

• Unsymmetric membrane-bending coupling 
stiffness terms and thermal coefficients 


Fig. 2 Traditional formulation of stiffened panels using 
beam and plate analogies. 
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STIFFENED PANEL EQUIVALENT PLATE 
FORMULATIONS 


Traditional Formulation Described 

Fig. 2 summarizes two analogies for formulating stiffened 
panel properties. The most simplistic analogy is to 
represent the stiffened cross section as detached beams. 
The neutral axis for the X-face is determined and the 
parallel-axis theorem is used to calculate bending stiffness 
on the X-face. The same process is performed on the Y- 
face. Coupling of forces and strains between the longitudi- 
nal and transverse directions is ignored by representing the 
panel as a series of separate beams. Therefore, the A 12 , 
D 12 , & (1-V 2 ) terms become omitted from the formulation. 

A more advanced approach is to represent the stiffened 
panel using plate formulation, Fig 2. This analogy 
correctly includes longitudinal-transverse directional 
coupling, but is still in error by representing the Y neutral 
axis as lying on the same plane as the X neutral axis, thus, 
omitting the unsymmetric, membrane-bending coupling 
response of the panel. 

The New Formulation Described 

Description of the new panel formulation begins with a 
review of the previously reported laminate formulations. 
For a layered material, the membrane, membrane-bending 
coupling, and bending stiffnesses are noted as 


(A A D ,) 


h/2 


/ f Q ij) k ( !, -z, z 2 )dz 

-h/2 


( 1 ) 


A:f, B:f, D,' 


h/2 


/ ( Q ij} k ( 1 ^ Z > z2 ) dz 


( 6 ) 


and 


h/2 


( Aj P<l , Bj P<l , Dj P<l ) = / (<E>:) k (1,-z, z 2 )dz (7) 


-h/2 


for the panel membrane, membrane-bending coupling, and 
bending stiffness terms and thermal coefficients. The 
asterisks indicate laminate and not lamina properties. The 
laminate properties are defined as 


Q 


>j 



( 8 ) 


$: = (9) 

1 t 

with material properties interpolated from a database, 
providing non-linear temperature and load dependent data 
based on the aircraft trajectory event. The FEA computed 
tension or compressive load and in-plane and through-the- 
thickness temperature gradients are used to generate these 
laminate or metallic sheet properties. If the panel sheets 
are laminates, the properties of the sheet are treated as 
being homogeneous, which is a valid assumption because 
the panel depth is much greater than the laminate 
thicknesses t x , t 2 , or t 3 . This assumption is also made in 
classical lamination theory because each ply is treated as 
being homogeneous even though it is a mixture of fiber 
and matrix. 


and reference 3 shows the corresponding thermal force 
and moment coefficients to be 


h/2 


( A", B", D") = / (0,) k (l,-z,z-)dz (2) 

In equation (1) the 


-h/2 


Qij - U:: T ^ 0) 


and are the transformed reduced layered elasticities of the 
laminae. [T] 4 is a fourth order transfer tensor and Q- are, 
as an example 


Q 


n 



(4) 


In equation (2) the 


O, = Q::«: T ^ (5) 


and are the transformed reduced layered thermal force 
coefficients, or are the material expansion coefficients. 
This approach extended to panel concepts has been shown 
in references 2, 3, and 6 to be 


The equivalent plate formulation of any stiffened panel 
shape, through extension of classical lamination theory, is 
accomplished by locating a reference plane, identifying its 
layers with a k, value, and defining the h heights from the 



Fig. 3 Laminate formulation is extended to stiffened 
panels. 
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reference plane. The panel layers, in this sense, are the 
facesheet and coresheet laminates and joining nodes. This 
approach produces the following general equations for 
panel stiffness terms and thermal coefficients 


A ; , -2B ; , 3D 


p\ - 


= E 

k=l 


/ u III 
( ^k-1 " 


u m \ 

h k ) - 

K / //''x 


( 10 ) 


(Qii) k 



w,. or 

k sin 9 


(m=l, m=2, m=3) 


Fig. 4 Any reference plane may be used to calculate 
properties by including unsymmetric data. 


( 11 ) 


/ 1 111 _ i UK 

{AT, - 2 B r, 3Dr ) = £ 1 k ' 1 " kj ($i )k 

k=l 3 


sin 9 


w k or — — - | (m=l, m=2, m=3) 


where S is the distance of the repeating pattern of 
corrugation and w is the width, t is the thickness, and 0 
is the angle of a stiffening segment (0=90° for perpen- 
dicular stiffeners). Each stiffness term and thermal 
coefficient is the summation of all laminate/metallic- 
sheet segments. If a segment is horizontal (in the plane 
of the panel) then its width is used. If a segment is 
vertical or at some angle, then its thickness divided by 
the sine of its angle is used. In this way, each segment 
and its shape can be accounted for in any panel concept. 

In equations (1, 2, 6, 7, 10, and 11) the [B] and {B a } 
are negative to match the sign convention of Fig. 3. 
These unsymmetric data are much more significant for 
stiffened panels than for laminates. A measure of the 
panel’s membrane -bending coupling can be visualized in 
Fig. 2. The distance between the X-face and Y-face 
neutral axes can only be accommodated by the [B] and 
{B“} data. Since the [B] and {B“} data permit a 
common X-face and Y-face reference plane, its location 
becomes a matter of choice. Conventional lamination 
theory uses the midplane, however, the vehicle outer 
mold line (OMF) surface, or an offset to account for a 
thermal protection system. Fig. 4, would be a more 
convenient choice. 

Kirchoff’s plane-sections-remain-plane hypothesis 

Fimitations and assumptions of the method as applied to 
any panel concept are the same as those applied in 
classical lamination theory. A primary assumption is 
that strain variation through the panel cross section 
follows the Kirchoff hypothesis for laminated plates. 
This hypothesis maintains that a normal to the midplane 
remains straight and normal upon panel deformation and 
that stresses in the XY plane govern the laminate 


behavior. Implications of this hypothesis are: 1) mem- 
brane strains vary linearly through the panel cross section 
2) stresses vary in a discontinuous manner through the 
cross section 3) the facesheet laminates are perfectly 
bonded to the coresheets, and 4) the bonds are infinitesi- 
mally thin and non-shear deformable. This implies that 
e p z , , y p xz , & Y yz = 0, in addition to the usual plane stress 
assumption of o p z , x p xz , & T p yz = 0. 

3-D FEAs (hat shape modeled with shell elements) were 
performed for the purpose of investigating the strain 
profile of the stiffened panel. Strain and curvature for 
thermal expansion and bending, caused by both in-plane 
(AT) and through-the-thickness (AG) temperature gradients 
are shown in Fig. 5. As shown with the bold line, the free 
edge of the FEM indicates a linear strain distribution 
through the panel depth. FEA of general mechanical loads 
shows the same linear strain profile. However, some 
loadings can cause high shear strain gradients of the webs 
near the panel boundary if they are free (with out joining 
structure). Even for these special cases, the non-linear 
strain profile was seen to extend only a short distance 
from the free edge. Since stiffened panels do not have 
large deflections relative to their depths, they do not 
appear to have the same geometrically non-linear concerns 
of "thin" plates. 



4036 clemen Is (75x) 


Fig. 5 FEA computed thermal strain profile through the 
panel depth. 


35th SDM Conference, AIAA-94-1579 


4 



Thermal load dependent residual strains and stresses 

Residual panel strains and stresses caused by thermal 
growth may be resolved in the same manner for varying 
stiffened shapes. They develop in stiffened panels, for 
each thermal loadcase, when the panel laminates elongate 
non-uniformly due to heating. The panel laminates cannot 
act independently and develop residual strains and stresses 
when forced to strain together as a unit of the panel. 
Panel curvature dictates that all laminate strains follow its 
through-the-depth strain profile. This profile is linear due 
to Kirchoff’ s hypothesis as discussed above. The residual 
strain is the difference between the strain that occurs in 
the laminate when made a segment of the stiffened panel's 
linear strain profile, and the strain that occurs in the 
laminate when allowed to thermally grow unattached to 
the panel. 

2-D FEA is able to use smeared equivalent plate 
properties for a stiffened panel because of this "plane - 
sections-remain-plane" hypothesis. In principle, during 2- 
D FEA, panel laminates strain together as a unit and thus 
do not include residuals. In order to quantify a laminate’s 
"design-to" strains, its residual strains must also be 
quantified and added to 2-D FEA computed strains, see 
ref. 3. A comparison of residual stresses from discretely 
meshed, 3-D FEA and residual stresses computed with 
coarsely meshed 2-D FEA that includes the method of 
computing residuals is shown in Fig. 6 (back page) 

Formulations Specific to a Panel Concept 

The application of the general stiffness formulation to a 
specific panel concept is presented below. Though the 
formulation is general and applies to any panel concept, 
some panel concepts have particular behaviors which must 
be addressed in unique ways. Special stiffness terms are 
identified for each unique behavior of a corrugated 
stiffened panel and included in the general equations (10) 
and (11), where required. 

The equations are grouped by panel direction in order to 
make them as brief as possible. A 13 , A-, 3 , B 13 , B-, 3 , Dj 3 , 
and D-, 3 are zero because the direction of corrugation is 
parallel to the panel's reference axis and because the 
laminates are balanced. 


X-axis; stiffnesses parallel to the corrugation 

The longitudinal stiffness terms for a corrugated stiffened 
panel are 



-2B 


p 

if 



(Qn )r (ho m -hD 


( 12 ) 


(Q u ) 2 , (h; n -hj n )Nw t 

S x 


2(Q U ) 2 (h‘ n -h ‘ n )t 2 2(Q U ) 2 (h‘ n -h 5 m )t 


sin0 S„ 


sin0 S„ 


where geometric variables such as coresheet angle 0, 
corrugation spacing S x , and widths of the coresheet top 
and bottom joining nodes Nw, and Nw b are shown in Fig. 
3. The laminate’s X axis is parallel to the panel’s X axis 
which is the corrugation direction. By using the OMF 
surface as the panel’s reference plane, the h variables are 
calculated as: h 0 =0, h^-tj, h 2 =h,-Nt t , h 3 =h r Nt/2, h 4 =-H/2, 
h 8 =-H, h 7 =h 8 +t 3 , h ( ,=h 7 +Nt b , and h 3 =h 7 +Nt b /2. The 
variables t lf t 2 , and t 3 are the top facesheet, coresheet, and 
bottom facesheet thicknesses. Nt, and Nt b are the thick- 
nesses of the coresheet top and bottom joining nodes. 
Higher panel bending stiffnesses will be calculated if the 
OMF is used but they will be balanced out with higher 
membrane-bending coupling stiffnesses. Note that a 
symmetric panel will now have non-zero membrane- 
bending coupling data. The advantages of using the OMF 
instead of the panel neutral axis are explained in reference 
3. 

The subscripts 1, 2, and 3 on the Cf terms represent the 
different isotropic materials or composite layups. 
Properties of these materials or layups are based on their 
non-linear temperature and load dependent data. The 
terms ( Cl ) 2 t and ( Q1 ) 2 b distinguish the coresheet top and 
bottom node laminate reduced stiffnesses from the middle 
coresheet reduced stiffnesses. Even though they are the 
same material or layup, a through-the-thickness tempera- 
ture gradient causes their properties to be dissimilar. 

Unlike the facesheets, the corrugated coresheet does not 
behave as a plate. Its nodes and mid portion strain in the 
longitudinal direction like a thin strip of plate or a beam. 
Because of this, the coresheet does not contribute to the 
panel stiffness terms 22 and 12 of equations (15) and (17). 
More subtle is the fact that the coresheet’s contribution to 
panel longitudinal stiffness as included in equation (12) is 
not effected by plate coupling as identified by equation 
(4). The coresheet’s accordion shape allows it to strain 
unconfined in the transverse direction eliminating the need 
for the familiar plate term (1-v 2 ). Coresheet laminate 
membrane stiffness (A u ) 2 which includes orthotropic 
coupling cannot be divided by its thickness to obtain 
( Q u ) 2 as shown in equation (8). Instead the coresheet’s 
uncoupled (Q u ) 2 terms equal the effective laminate 
engineering elasticity, E x 


E.. - 


A ll( 1 -'Vyx) 


or more generally for unsymmetric 
laminates 


E 


X 


1 

An t 


(13) 


or unbalanced 


(14) 


(Qn) 2 - (h6 n -h 7 m )Nw b 

s x 


(m=l, m=2, m=3) 
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Y-axis; stiffnesses perpendicular to the corrugation 

The transverse stiffness terms for a corrugated stiffened 
panel are 

(15) 

(A'’, -2B 2 2 , 3D'’) = 


(Q22 )t (ho n -hi n V 


(Q 22 (hr-h 2 m )Nw t 

s x 


(m=l, m=2, m=3) 


As noted in the longitudinal stiffness equations, the 
coresheet does not contribute to the transverse stiffness. 
However, if the coresheet top node is securely attached 
to the facesheet, as in the case when it is relatively wide 
and brazed to it, then its term shown in equation (15) is 
significant. If so, then its contribution to the longitudi- 
nal stiffness as shown in equation (12) is based on 
equation (8) and not equation (13) or (14). If the top 
coresheet node is spot welded to the facesheet then it 
should be omitted from equation (15). The transverse 
stiffness then becomes dependent on the fabrication 
method and must be properly considered with either 2-D 
or 3-D FEM’s. 

Curvature of the corrugated stiffened panel in the 
transverse direction causes coresheet flexure as noted in 
Fig. 7 As depicted, pure bending caused by an applied 
rotation on the end produces the same rotation at the 
facesheet and coresheet joint. This rotation causes the 
hat to take on the deflected shape shown, which is 
resisted by both coresheet and facesheet flexure. This 
additional D)’, bending stiffness can be quantified by 
taking ratios of the coresheet flexure stiffness to the 
facesheet flexure stiffness. With new shape variables b 
and b 2 . Fig. 8, the additional stiffness term is 




Panel numbing 
cauc coresH e f t 
flexme 


Fig. 7 Transverse panel bending is resisted by both the 
facesheet and coresheet. 



Fig. 8 Detailed shape variables are required for 
formulations. 



-(^22)2 b + (^22)2'’ b (D 22 )ib 

(D 2 2)ib 2 + (D^NwJ (S x ) 


which gives the panel over 180% more stiffness for the 
example panel shape shown later when the facesheet 
and coresheet are the same layup. In this equation, the 
laminate’s bending stiffnesses are used instead of the 
Q7 (h m k , - h m k ) term to obtain more accuracy. The 
laminate’s bending and membrane-bending stiffnesses 
could be used in equations (15) and (17) as well. 

In summary, equation (15) is dependent on the choice of 
reference plane. It includes both parts of the parallel 
axis theorem which is the rotation and extension of the 
facesheet in the cross corrugation direction. Equation 
(16) is not effected by choice of reference plane, as it 
only quantifies the coresheet resistance to rotation. 


XY-axis; directional coupling stiffnesses 

The directional coupling stiffness terms for a corrugated 
stiffened panel are 



-2B 


p 

12 ’ 



(17) 


(Q12 )t (ho n -h;V 


(Q 12 > 2 ' (hr-h 2 m )Nw t 
S x 


(m=l, m=2, m=3) 


As noted in the transverse stiffness equations, if the 
coresheet top node is securely attached to the facesheet, 
then its term shown in equation (17) is significant. If 
the top coresheet node is not secure to the facesheet 
then it should be omitted from equation (17). 
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XY-axis; in-plane shear, bending-twisting, and their 
coupling stiffnesses 

The hat-shaped, closed cells of the corrugated stiffened 
panel behave like torque tubes. As such, they provide 
considerable additional stiffness to the panel B 13 ^, 

DP33 terms. These terms must be accurately formulated 
in order to take advantage of the unique characteristics of 
the corrugated panel. The A p , , and D p 3 3 formulations 
are variations of those found in references 8 and 9. The 
B P 33 formulation is new. The D p 33 formulation has been 
improved. 

Formulation of these 33 terms is dependent on the way the 
panels will be joined together. A typical fuselage panel 
and ringframe assembly is shown in Fig. 9. This design 
fully attaches the panel facesheets to the ringframe with a 
splice plate and many fasteners. Therefore, the three 
translational and the three rotational facesheet boundary 
conditions are linked together. The coresheet is connected 
to the ringframe flange with a fastener as depicted in the 
cut away view. This attachment scheme allows the out-of- 
plane rotational degree-of-freedom of the bottom coresheet 
to be unconstrained, permitting limited hat cell waipage. 
Listed below are stiffness terms for this common design. 
Once the 33 terms have been calculated based on their 
application, they can be entered in the 2-D FEM without 
having to manually release this rotational degree-of- 
freedom as required for a 3-D FEM. Fig. 10 shows the 
proper panel shear deformation and Fig. 11 the proper 
bending-twisting curvature, both with hat warping, as 
performed with 3-D FEA. 


closed cell. Quantifying the ratio (q cr ) of the cell shear 
flow, to the panel shear (N ) aids the A p , , formulation. 
Calculation of this ratio is indeterminate but can be solved 
by using the panel strain kinematics and the cell shape 
variables 


(19) 




(A33) 1 


2b, 


Nw,, 


(^33)2 (^33)2 


where, for brevity, the following variables are used 


(A33) 1 
(^33)2 
(^33)2'’ 

0^33)2' 

(•^33)4 

and, from Fig. 8 
b = 

b 2 = 


( 20 ) 

(0.33)1 (V h i) 

(Q 33 )2 (2 
(Q33 )2 b (^6 _ ^7) 

(0.33)2' (^ 1 _ ^ 2 ) 

(033)1 (V h i) + (0.33)2' (V h 2 ) 

(21) 

S x - Nw t 

h 3 - h S 

sin0 



Fig. 9 A typical panel joint design assembly is used for 
formulations. 


The in-plane shear stiffness term for a corrugated 
stiffened panel is 

(18) 


N w t + 2 b, q cr + Nw b q cr 
(^33)4 (^33)2 (A. 33 ) 2 b 


Shear flow and laminate shears are 

q = 

K ~ N x>, - N xy, + q 


( 22 ) 



Fig. 10 Computed shear deformation of a 3-D FEM. 


This equation is derived from the knowledge that a 

constant shear flow (q) exists around the hat-shaped. Shear flow around the cell causes coupling between in- 

plane shear and bending-twisting. This response is 
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included in B 1 ’ , ,. BP 3 , has two parts: the first part noted 
as B quantifies the twist caused from the cell shear 
flow, and the second part quantifies the panel twist from 
the in-plane shear of the facesheet as a function of the 
reference plane location 

(23) 


B p 


1 In ' x n, 2 h h (Q ' ' >2' ( h i 2 - h 2 2 ) Nw : 
'7 (Q.y )i ( h o ) + - 


+ BiV 


where 


B .5 = a q ^ 


and the cell area and height is 


A = 


Nw, 


H 


cell 


H 


cell 


h n - — - h. 


(24) 


(25) 


The derivation of B [, begins by noting that the cell 
shear flow causes torque in proportion to the area 
enclosed by the cell. The Bredt formula quantifies this 
as 


T = 2 A q 


(26) 


The relation of torque to the bending-twisting moment is 


M 


\y 



(27) 


The cell shear flow can be written in terms of 


q = Ter A .y 

by knowing that 

N p - A p v p 

iN xy - /VyYxy 


(28) 

(29) 


and using equation ( 22 ). 
shear flow can be written as 



The coupling term due to 


(30) 

M p 

xy 


when all other strains and curvatures are zero. By 
substituting equations (26), (28), and (29) into equation 
(27), equation (30) becomes equation (24). 

Applied moment M xy or torque is resisted by the panel's 
bending-twisting stiffness D 1 ’, ,. The corrugated panel's 
closed cells provide much more rigidity than an open 
cross section such as the blade stiffened panel. The 
resulting differences are important because the D 33 term 
has four times the effect of increasing the panel’s 
buckling load than does the primary bending stiffness 

D„. 


DP 33 also has two parts: the first part noted as D [, 

quantifies the resistance to twist caused from the cell shear 
flow, and the second part quantifies the panel resistance to 
twist from the in-plane shear of the facesheet as a function 
of the reference plane location 

(31) 


rx p 1 [,^ - .1 ^ ,q (O.v, ) 2 ' (hi’-h/) Nw : 

= - (Qy \ (ho -hi ) + 


n p 

-Gy 


where 



A 2 



(32) 


is written in terms of M and ds is the length of the 
cell perimeter. The D [, formulation is a form of the 
familiar closed cell torsional rigidity equation. It is 
applicable to any closed shape and is independent of the 
reference plane location. Applied to the hat cell, 
equation (32) becomes 



G f + 


AA 

A 


2 b 2 
(Ay) 2 


+ 


Nw b 

(Ay) 2 b 


(33) 


where G f accounts for the cell facesheet shear flexibil- 
ity. G f -b/ (A„), when the cell is not part of a panel 
and the entire cross section is allowed to warp. G f = 0 
for the case where the cell is still not part of a panel, 
but the facesheet is prevented from warping, as in the 
design of Fig. 9. When the cell is part of a panel, the 
shear flexibility of the facesheet between cells allows 
shear flow into the cell facesheet even though the 
facesheet’ s average warpage and shear strain is pre- 
vented. Therefore, for the hat panel, 

G f - Nw t /S x (Ay) 4 



Wole: 



Fig. 11 Computed bending-twisting curvature of a 3-D 
FEM. 
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STIFFENED PANEL EQUIVALENT PLATE 
VERIFICATION 


Twelve unique stiffness terms 

An anisotropic panel has eighteen unique stiffness terms. 
The 13 and 23 terms of the [A], [B], and [D] matrices 
become zero leaving twelve unique terms for a corrugated 
panel that is fabricated with balanced layups and is 
modeled with the panel’s axis parallel or perpendicular to 
the corrugation. Values for these twelve unique stiffness 
terms were calculated using the formulations of this paper 
and then compared to 3-D FEA. Two reference planes 
were used to test the generality of the equations. Table 1 
lists values for the panel midplane and Table 2 lists values 
for the facesheet midplane. The facesheet midplane is a 
useful reference plane because it is close to the OML 
surface and because it is the location where the smallest 
values for B p 12 , B p ~, 2 , D P j 2 , and D p 77 occur. 

The twelve stiffness terms were determined independently 
of each other by applying one known strain or curvature 
and setting all others to zero. By isolating each panel 
response this way, the six equations and six unknowns that 
relate panel strain and curvature to load 


l N 'l 

A 

B s' 

ht 

1 M l 

B 

D 

hi 

l 1 J 6x1 

L ‘J 

y_ 

6x6 l 1 J 


were used for measuring each stiffness term with 3-D 
FEA. 

Six cases per reference plane (three unit strains and three 
unit curvatures) were applied to two different 3-D FEMs: 
a large panel area with many hats and a small panel area 
with many elements to represent the hat shape. Unit 
strains and curvatures were applied with MSC/- 
NASTRAN™ RBE2 rigid elements on each side of the 
panel. Unit enforced displacements or rotations were 
applied one at a time, maintaining all other deformations 
to be zero. The resulting reacted forces and moments 
were then used to measure any applicable terms. For 
instance, for a unit e x , four terms can be measured 

A u =N x , A 12 = N y , B u = M x , B 12 = M y (35) 

The 3-D FEMs used for these comparisons had an 
inherent shortcoming of having the shell elements of the 
top coresheet co-planar to the facesheet elements. Co- 
planar elements were necessary to maintain connectivity 
to the same grid points. In lieu of off-setting the top 
coresheet properties in the FEM, the hj distances shown in 
Fig. 3 were adjusted for these comparisons. Also, for 
these comparisons, the stiffness terms were calculated 
assuming the panel was fabricated in a way that did not 
provide secure enough attachment between the top 
coresheet and facesheet to allow the transverse properties 
of the top coresheet to be used in any terms. 


Panel Design and Temperatures 

The panel cross section shape, dimensions, and laminate 
layups used for the stiffness term comparisons is shown in 
Fig. 12. To the left of the section is a temperature profile 
which is typical of those analyzed for hypersonic flight. 
The shaded rectangle represents a uniform in-plane gradi- 
ent of 555°F (625°F-70°F). The double shaded triangles 
represent a through-the-thickness gradient of 300°F/in. By 
superimposing the two gradients, the variation of temper- 
ature through the panel’s depth is known, as illustrated 
with the bold line. The facesheet’s average temperature 
of 842.5°F. and the coresheet’s average temperatures of 
832.75°F„ 617. 5°F., and 402.25°F. were used for 
interpolating the material database. Temperature 
dependent laminate data Q7 and were used for 
formulating panel stiffness terms and thermal coefficients, 
in order to make the comparisons more challenging. The 
temperature gradients were not applied.as loads, because 
their induced thermal forces would have made the 
solutions indeterminate. 



Room temperature = 70° F Metal matrix composite material 

In-plane temperature gradient (AT= 555° F) 

Through-the-thickness temperature gradient (AG = 300° F/in) 

l/FCSHT2/Collier 


Fig. 12 The corrugated stiffened panel and temperature 
profile used for analysis. 


Observations 

Panel in-plane shear, bending-twisting, and their couplings 
are highly dependent upon boundary conditions, thus, 
equivalent plate formulations and 3D FEMs must be 
tailored appropriately. The A p ,^, B p and D p ,, are 
greatly increased by attaching the hats to a stiff support 
such as a ringframe. The A p ,^, B p ,^, and D p , , formu- 
lated values based on the design in Fig. 9 matched within 
0.1% of FEA performed with the finely meshed model. 
The remaining nine stiffness terms were computed only 
with the larger, coarser meshed 3D FEM. They too 
compared well with formulated values (within 1%) and are 
listed in tables (1) and (2). In short, all twelve stiffness 
terms are correctly formulated for the hat-shaped, 
corrugated stiffened panel. 

The close comparisons to measured values indicate that 
this particular equivalent plate approach provides the 
capability to match 3-D FEA computed response for any 
general thermomechanical load condition. Implication of 
this match is that proper equivalent plate formulations of 
stiffness and thermal coefficients provide the capability to 
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quantify the thermoelastic response of aircraft using 
relatively coarse meshed 2-D models. 

Regardless of choice of reference plane, the same correct 
solution of laminate strains and stresses are obtainable, 
even though the calculated stiffnesses and thermal 
coefficients, and computed moments vary greatly. 
Reference plane dependent, FEA computed 
thermomechanical N P j and M P j are multiplied by the 6x6 
inverted stiffness matrix, see equation (34), to resolve 
reference plane e P j and k P j. These panel reference plane 
8 P j and k P j are then used in equation (36) to quantify 
strains at any location in the panel depth. A laminate 
strain is the summation of the panel's reference plane 
strain and the additive or subtractive contribution of the 
curvature 

<«h} k = e p j - kF H k (36) 

where H k is the distance between the reference plane and 
the panel depth location k. 


CONCLUSIONS 

Formulations presented in this paper provide the capability 
to model stiffened composite panels of any cross sectional 
shape with a single plane of shell finite elements. Twelve 
unique stiffness terms were formulated for a corrugated, 
hat-shaped composite stiffened panel by extending 
classical lamination theory to the stiffened cross section, 
and by including the unique characteristics of the closed 
cells. A new formulation for the bending-twisting 
coupling with in-plane shear stiffness term (B p 33 ) was 
derived and an improved D p 33 was presented. The 
formulated values were compared with 3-D FEA. All 
terms had differences less than 1%. Implication of this 
agreement is that the equivalent plate formulations of 
stiffness and thermal coefficients presented in this paper 
provide the capability to quantify the thermoelastic 
response of aircraft using relatively coarse meshed 2-D 
models. 
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Table 1. Comparisons of formulated stiffnesses to 3D FEA computed stiffnesses, 
(the reference plane is the panel midplane) 


CASE 

A u 

A)2 

A 22 

A 33 

B u 

B)2 

B 22 

B 33 

D n 

D)2 

D 22 

B^33 

FORMU- 

14.1 

2.60 

6.97 

2.12 

- 5.71 

- 1.91 

- 5.12 

- 1.34 

6.14 

1.40 

3.75 

1.24 

LATED 

(+05) 

(+05) 

(+05) 

(+05) 

(+05) 

(+05) 

(+05) 

(+05) 

(+05) 

(+05) 

(+05) 

(+05) 

3D FEA 














14.1 

2.60 



- 5.70 

- 1.91 







8 y 


2.60 

6.97 



- 1.91 

- 5.11 






Yxy 




- 




- 





K x 





- 5.73 

- 1.91 



6.15 

1.40 



K y 






- 1.91 

- 5.12 



1.40 

3.75 


K v „ 








. 




. 

xy 














Table 2. Comparisons of formulated stiffnesses to 3D FEA computed stiffnesses, 
(the reference plane is the face sheet midplane) 


CASE 

A u 

A)2 

A 22 

A33 

Bn 

B 12 

B22 

B33 

D„ 

D)2 

D 22 

B^33 

FORMU- 

14.1 

2.60 

6.97 

2.12 

4.63 

0.00 

0.02 

1.28 

5.35 

5.41 

1.54 

1.59 

LATED 

(+05) 

(+05) 

(+05) 

(+05) 

(+05) 

(00) 

(00) 

(+04) 

(+05) 

(+01) 

(+02) 

(+04) 

3D FEA 














14.1 

2.60 



4.63 

0.00 







e v 


2.60 

6.97 



0.02 

0.00 






Yxy 




2.12 




1.28 





Kx 





4.65 

0.00 



5.38 

5.53 



K y 






0.00 

0.02 



5.43 

1.54 


K 








1.29 




1.59 

xy 














3000 




Land 


0 

-1000 



psi 




Fig. 6 Comparison of theoritical thermal residual stress to a 3D FEA solution. 
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